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Abstract 

It has been suggested in Ref. [A. E. Bondar and V. L. Chernyak, Phys. Lett. B 612, 215 (2005)] 
that the disagreement between theoretical calculations and experimental observations for the rate 
for the process e + e~ — > J/ifj+r/ c at the B factories might be resolved by using the light-cone method 
to take into account the relative momentum of the heavy-quark and antiquark in the quarkonia. 
The light-cone result for the production cross section in Ref. [A. E. Bondar and V. L. Chernyak, 
Phys. Lett. B 612, 215 (2005)] is almost an order of magnitude larger than existing NRQCD 
factorization results. We investigate this apparent theoretical discrepancy. We compute light-cone 
distribution functions by making use of quarkonium wave functions from the Cornell potential 
model. Our light-cone distribution functions are similar in shape to those of Ref. [A. E. Bondar 
and V. L. Chernyak, Phys. Lett. B 612, 215 (2005)] and yield a similar cross section. However, 
when we subtract parts of the light-cone distribution functions that correspond to corrections of 
relative-order a s in the NRQCD approach, we find that the cross section decreases by about a 
factor of three. When we set certain renormalization factors Zi in the light-cone calculation equal 
to unity, we find a further reduction in the cross section of about a factor of two. The resulting light- 
cone cross section is similar in magnitude to the NRQCD factorization cross sections and shows 
only a modest enhancement over the light-cone cross section in which the relative momentum of 
the heavy-quark and antiquark is neglected. 
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I. INTRODUCTION 



One of the outstanding issues in quantum chromodynamics (QCD) in recent years is the 
discrepancy between the theoretical predictions and the experimental measurements for the 
exclusive production of J/ip + r\ c in e + e~ collisions at the B factories. The cross section 
times the branching ratio into at least two charged tracks has been measured by the Belle 
collaboration to be 25.6 ± 2.8 ± 3.4 fb (Ref. [1, 2]) and by the BaBar collaboration to be 
17.6 ± 2.8^2A fb (Ref. [3]). In contrast, calculations by Braaten and Lee [4] and by Liu, 
He, and Chao [5] using the nonrelativistic QCD (NRQCD) factorization approach [6] at 
leading order in a s yield predictions of 3.78 ± 1.26 fb and 5.5 fb, respectively. 1 A similar 
disagreement between NRQCD factorization and experiment holds for production of Xco 
and rj c (2S) mesons in conjunction with a J/ip meson. A recent calculation of corrections of 
next-to-leading order in a s leads to an enhancement of the theoretical prediction of about a 
factor 1.8 (Ref. [7]). Nevertheless, a substantial discrepancy between theory and experiment 
remains. 

It has been suggested in several recent papers [8-11] that predictions for the J/t/j + r\ c 
production cross section in the light-cone formalism might be in better agreement with 
experiment. In the present paper, we focus on the work of Bondar and Chernyak in Ref. [8], 
henceforth referred to as BC. In BC, it is claimed that an enhancement of the cross section 
occurs in the light-cone formalism because the expressions for the cross section contain 
an integration over the light-cone momentum. In contrast, in the NRQCD factorization 
approach, cross sections are computed as an expansion in the velocity v of the heavy quark or 
antiquark in the quarkonium rest frame. In the leading order in v, the relative momentum of 
the heavy quark and antiquark are neglected, which is equivalent in the light-cone formalism 
to approximating the light-cone distribution amplitude (light-cone wave function) by a delta 
function at a light-cone momentum fraction of 1/2. Effects of the nonzero relative momentum 
of the heavy quark and antiquark in NRQCD are taken into account through terms of higher 
order in the v expansion. 

The NRQCD factorization and light-cone formalisms are both believed to be valid approx- 
imations to QCD for the process e + e~ — > J/ip + r\ c in the limit in which the hard-scattering 

1 The differences in the predictions arise from different choices for a s and the charmonium wave functions 
at the origin and from the inclusion of QED effects in the calculation of Ref. [4] . 
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momentum transfer is much greater than either the QCD scale Aqcd or the charm-quark 
mass m c . 2 In NRQCD factorization, one makes the further assumption that the cross section 
can be written as an expansion in powers of v. Since the light-cone and NRQCD factoriza- 
tion approaches are both first-principles methods, their predictions should be in reasonable 
agreement unless, as is asserted in BC, the v expansion of NRQCD breaks down for the 
process e + e _ — > J/ip + i] c . 

In this paper we investigate several issues with regard to the calculation in BC by making 
use of light-cone distribution amplitudes that are derived from potential-model quarkonium 
wave functions. We first compare the light-cone distribution amplitude that we derive with 
the model light-cone distribution amplitudes that are used in BC. We find that they are 
quite similar in shape, but with significantly different asymptotic behavior at large light- 
cone momentum fraction, and that they lead to similar values for the cross section for 
e + e~ — > J ftp + r] c . We then argue that a part of the contribution to the cross section 
that is computed in BC comes from the high-momentum tail of the light-cone distribution 
amplitude. Such high-momentum contributions are taken into account in NRQCD through 
contributions to the cross section of relative order a s and higher. Furthermore, since our 
potential-model wave function is accurate only for momenta that are much less than m c , 
it is not legitimate for us to include the high-momentum contribution in our light-cone 
calculation. When we subtract this contribution, we find that the cross section is reduced 
by about a factor of three. A further reduction of about a factor of two occurs when we 
replace with unity certain renormalization factors Zi that appear in the calculation in BC, 
but which have no counterpart in a conventional NRQCD factorization calculation. With 
these modifications, the light-cone calculation gives a cross section that is comparable to 
that of the NRQCD factorization calculations. 

The remainder of this paper is organized as follows. In Sec. II, we describe the Cornell 
potential model that we use. In Sec. Ill, we present the relationships of the required light- 
cone distribution amplitudes to the Bethe-Salpeter (BS) wave function, use the BS equation 
to compute the BS wave function in terms of a potential-model wave function, and give the 
correspondences between the terms in the resulting BS wave function and the required light- 
cone distribution amplitudes. Sec. IV contains the detailed expressions for the light-cone 

2 For a discussion of the current status of proofs of NRQCD factorization, see Refs. [12-15]. 
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distribution amplitudes in terms of the potential-model wave function and a re-arrangement 
of these expressions into forms that are amenable to numerical calculation. In Sec. V, we 
identify the high-momentum part of the light-cone distribution amplitude that is properly 
treated as an order-a,, correction to the cross section. We also compute and compare the 
asymptotic forms of the light-cone distribution amplitude and the high-momentum part in 
the high-momentum limit. In Sec. VI, we present the expression for the cross section in the 
light-cone formalism, give our numerical results for the cross section at £>-factory energies, 
and discuss the implications of these results. Finally, in Sec. VII, we summarize our findings. 



II. POTENTIAL MODEL 



In this section, we describe the calculation of the potential-model wave function. We 
note that, if we knew the heavy-quark potential exactly, then we could calculate the heavy- 
quarkonium wave function in a potential model, up to corrections of relative order v 2 
(Ref. [16]). We make use of the Cornell potential model of Ref. [17]. The Cornell po- 
tential provides a reasonably good fit to heavy-quark potentials that are measured in lattice 
calculations. 3 

The Schrodinger equation for the radial wave function for a quark-antiquark (QQ) pair 
with orbital-angular- momentum quantum number £ in a central potential V(r) is 

R(r)=e B R(r), (1) 

where 6b is the binding energy, r is the distance between the heavy quark and antiquark, 
and m is mass of the heavy quark. Note that the reduced mass of the pair is \i = m/2. 
The Cornell potential [17] is given by 

V(r) = -- + - 2 , (2) 
r a z 

where the parameters k and a determine the strength of Coulomb and linear potentials, re- 
spectively. For a color-singlet QQ pair, the Coulomb-strength parameter k can be expressed 
in terms of an effective strong coupling a s as 

k = a s C F , (3) 

3 For a recent review that discusses heavy-quark potentials from lattice measurements, see Ref. [18]. 
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where Cf = 4/3. The parameter a is related to the string tension a as 



o = - 2 - (4) 



Following Ref. [17], we replace k and r by dimensionless variables A and p: 

K = (WT 2/3 A, (5a) 

r = p = fl(ma)" 1 '' 3 p. (5b) 
ran 

Substituting the Cornell potential into Eq. (1), we can rewrite the radial equation: 

u{p) = 0, (6) 



d 2 eu + i) x 

+--P+C 



dp 2 p 2 p 

where the dimensionless function u(p) and dimensionless energy eigenvalue £ in Eq. (6) are 
related to their physical counterparts as 

R(r) = (7a) 
V a z p 

e B = m{ma)~ 4/3 (. (7b) 
The wave functions are normalized according to 

/>oo />oo 

/ \u(p)\ 2 dp= / \R(r)\ 2 r 2 dr = 1. (8) 
Jo Jo 

In Eq. (2.19) of Ref. [17], the wave function at the origin is expressed in terms of the 
model parameters m and a as 

\m\ 2 = £^ (i + Mp~ 2 )) , 0) 

where the expectation value (p~ 2 ) is defined by 

/ ^ r Hp)} 2 , , x 

(p > = / U S ZLrf P- ( 10 ) 

Jo P 

Numerical values of (p~ 2 ) for various values of A are given in Table I of Ref. [17]. As we 
explain below, we fix the value of 1^(0)1 by using the leptonic decay width of the J/if). We 
need an additional constraint to fix both of the parameters m and a for a given A. Following 
Ref. [17], we use the mass splitting between the J/ip and the ^>(2S). Eq. (2.12) of Ref. [17] 
expresses the mass splitting as 

M^s) - Mj/^ = m(ma)- 4 / 3 (C 20 - Go), (11) 



6 



where C«o is the eigenvalue for the radial equation (6) for the principal quantum number n 
and the orbital- angular-momentum quantum number £ = 0. The values for £10 as a function 
of A are given in Table I of Ref. [17]. We use Eq. (A4) and Table II of Ref. [17] to determine 
C20 for various values of A: 

C20 = 4.0879 - 0.5826A - 0.0302A 2 . (12) 
Solving Eqs. (9) and (11), we can determine the model parameters m and a: 

m - C ™ ~ Cl ° ( 4 ^ )I 2 V /3 (13a) 
m ~ M+w-Mj/oKl + Xifi- 2 )) ' [ } 

From Eqs. (3) and (5a), it can be seen that the effective strong coupling a s for the bound- 
state potential as a function of A is given by 

«. = £~ (ma)" 2 / 3 . (14) 

As we have mentioned, we determine [^(O)! from the leptonic decay width of the J /if). 
NRQCD matrix elements for J/ip that correspond to [^(O)] 2 are given in Table I of Ref. [4]. 
These values are (0) s = 0.208 GeV 3 for LO and 0.335 GeV 3 for NLO, where the identifiers 
LO and NLO indicate that the short-distance coefficient for the leptonic decay rate of the 
J/ijj has been computed to leading order or next-to-leading order in the strong-coupling 
constant. Expressing [^(O)] in terms of these matrix elements, we obtain 



{ Ohl* = I 0-18619 GeV 3 / 2 (LO), 
2N c 1 0.23629 GeV 3/2 (NLO). 

In this paper, we use the LO value of [^(O)!, which corresponds to the formula 

r[j/^^e + n = ^^\m\ 2 , (is) 

m c 

where e c = 2/3 is the fractional electric charge of the charm quark. If one uses the NLO 
value of 1^(0)1, then the average momentum-squared of the wave function is reduced by 
about 10% (Ref. [19]). Therefore, we do not expect that the use of the NLO value would 
change our results qualitatively. 

Using |^(0)| = 0.18619 GeV 3 / 2 , M J/f = 3.096916 GeV, M f{2S) = 3.686093 GeV, Table I 
of Ref. [17], and Eq. (12), we can compute m and a from Eq. (13). The results for various 



TABLE I: Potential-model parameters and derived quantities as a function of the strength A of the 
Coulomb potential. The definitions of the parameters and derived quantities are given in the text. 
The parameters are computed using the inputs |Y>(0)| = 0.18619 GeV 3 / 2 , Mj^ = 3.096916 GeV, 
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values of A are shown in Table I, along with values for a s from Eq. (14), a from Eq. (4), and 
7c, which is related to the binding energy for a pure Coulomb potential and is defined by 

7c = ^ma s C F . (17) 



Lattice measurements of the heavy-quark potential yield values for effective coupling a s 
of 0.22 in the quenched case and approximately 0.26 in the unquenched case [18]. A lattice 
measurement of the string tension K = a (Ref. [20]) gives Ka\ = 0.0114(2) at a lattice 
coupling j3 = 6.5, where ol is the lattice spacing. Lattice calculations of the hadron spectrum 
at = 6.5 yield values for l/a L of 3.962(127) GeV (Refs. [21, 22]) and 3.811(59) GeV 
(Refs. [21, 23]). These yield values of the string tension of K = 0.1790 ± 0.0119 GeV 2 and 
K = 0.1656±0.0059 GeV 2 , respectively. Comparing the results of these lattice measurements 
with Table I, we conclude that A = 0.7 is a reasonable choice for the value of the potential- 
model parameter, and we use that value in our analysis. 

We express the radial equation (6) as a difference equation, and integrate it numerically. 
The result, for the choice of parameters in Table I that corresponds to A = 0.7 is shown in 
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p 

FIG. 1: The dimensionless wave function u(p) as a function of the dimensionless variable p. Input 
parameters are m = 1.14710 GeV, a = 2.38139 GeV" 1 , and A = 0.7. The wave function at the 
origin is taken to be |^(0)| = 0.18619 GeV 3 / 2 , which is the value for the J/tp wave function at the 
origin that is designated as LO in Ref. [4]. 

Fig. 1. 

III. RELATIONSHIPS OF THE LIGHT-CONE DISTRIBUTION AMPLITUDES 
TO THE POTENTIAL-MODEL WAVE FUNCTION 

A. Light-cone distribution amplitudes 

In this subsection, we re-produce the expressions that are given in BC for the light-cone 
distribution amplitudes in terms of the BS wave functions. 

For the J/ip with polarization e, the light-cone distribution amplitudes are defined by 

+ Z^*)^ ^(^) + /r^je^ ^ W ^a(^) ) • (18) 

M J/* J a/3 
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Here Q and Q are the heavy quark and antiquark fields, respectively, fi* is the renormal- 
ization scale, a^ u = | [7^,7^], and z± — z and z 2 = 1 — z are the fractions of the meson 
momentum P + = P° + P 3 = go carried by the quark and antiquark, respectively, at large 
P 3 . In deriving these expressions, the authors of BC neglect the quark transverse momen- 
tum inside the heavy quarkonium in comparison with the quark mass. This is a further 
approximation that goes beyond the standard light-cone formalism. As is pointed out in 
BC, the light-cone distribution amplitude V(z) may be eliminated in favor of the standard 
leading- twist wave function V L (z) of a longitudinally polarized J/t/j (Ref. [24]): 

(J/i>x=o(P)\Q(xhM-x)\U) = fvP, f d Zl e w < z ^V L {z). (19) 

Jo 

Taking e to correspond to helicity A = in Eq. (18) and using Mj^e^ =0 — > P M for large P, 
one obtains 

V L (z)=V ± (z) + V(z). (20) 

For the ?7 c -meson, one has 

(Vc(P)\Q(3(x)Q a (-x)\0)„* = j\z^< n -^^-P A {z)-^) lb P P {z) 

+fl{^)a, v P^x v P T {z)\ . (21) 

) a/3 

B. Relationships to potential- mo del wave functions 

In order to express the light-cone distribution amplitude in terms of a potential-model 
wave function, we now derive an approximate form for the BS wave function in terms of a 
potential-model wave function. Here we make explicit use of an expansion in the velocity v 
of the Q or Q in the quarkonium rest frame. Of course, such an expansion is not valid at 
large QQ relative momenta. However, as we will discuss in detail in Sec. V, we ultimately 
subtract such large-momentum contributions in the cross-section calculation. 

Consider a heavy QQ pair bound in a potential V. The four-momenta for the pair can 
be written in terms of quarkonium momentum P and their relative momentum p as 

PQ = l P + P, ( 22a ) 
PQ = \P ~ P- (22b) 
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In the rest frame of the heavy quarkonium, P and p are given by 



P = (M meson ,0) = (2m + e B ,0), 

V = ip°,p), 



(23a) 
(23b) 



where e-Q is the binding energy of the quarkonium. Let ip(x) be the potential-model wave 
function satisfying the time-independent Schrodinger equation: 



(24) 



where \i = to/2 is the reduced mass and p is the momentum operator. We specialize to the 
S-wave case. Then ip(x) — ip(r), where r = \x\. 

We can form a trial BS wave function, accurate to leading order in v, by appending spin 
wave functions onto ip. In the rest frame of the quarkonium, the spin-singlet and spin-triplet 
wave functions are 



n U = £%(%)« = ±= [(l + 7°)7 5 ] a6 = ^ 
i,j V V * 

iitu = = ~ [(i+7°) i] ab = ~ 



1 + — 

2to 



(25a) 
(25b) 



ab 



where £ and fj are Dirac spinors for a static free quark and antiquark, respectively. Cij( l So) 
and Cij( 3 Si) are Clebsch-Gordan coefficients of the spin-singlet and spin-triplet states for the 
spin states i and j of Q and Q, respectively, and a and b are the Dirac indices. SU(3) color 
indices are suppressed. In the last equalities of Eqs. (25), we have written the expressions 
in a Lorentz-invariant form. 4 Our trial wave function is then ■0(r)(n A ) ab , where A takes on 
the values or 1. In momentum space, the trial BS wave function is 



iP(p)2n5(p°)(Tl x ) ab = J d 4 xe-^(r)(U x ) ab . 



(26) 



The momentum-space form of the BS equation [28] is 
d 4 q 



*dT(p) 



(27^ 



S?(p) \K{p- q )] a , a ,, b „ b , ^ BS (g) ,„„ S%) 



J a"b" 



J b'b 



(27) 



4 These Lorentz-invariant expressions were first given in Refs. [25, 26]. Exact expressions, correct to all 
orders in v were given in Ref. [27]. 
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where Sp is the heavy-quark propagator. In our case, we can take the interaction kernel K 
to be 

[K(P - q)] a . a ..,w = la>a»HV(p - qH, v , (28) 

where 

V(p-q) = J d 3 xe^ p - q > x V(r) (29) 
is the Fourier transform of the nonrelativistic potential. Then we can rewrite the BS equation 



as 



*5?(p) 



(2vr) 4 p 2 Q - m 2 + ie 



(30) 



Since we have neglected retardation effects, which are of relative order v 2 , the interaction 
kernel (28) is independent of (p— g)°. Consequently, the integration over g° on the right-hand 
side of Eq. (30) produces the equal-time BS wave function. Integrating Eq. (30) over p° to 
form the equal-time wave function on the left-hand side, we obtain the Salpeter equation [29]: 



where 



dp f d 3 q [(^Q + m) 7 ] aa , 



2n J (27r) 3 Pq — m 2 + ie 



*L(p) 



iV(p-q) * s a , v (q 



m 2 + ie 



(31) 



(32) 



is the equal-time BS wave function, which is also known as the Salpeter wave function. 

We can obtain a wave function that is accurate through corrections of order v by substi- 
tuting the trial wave function (26) into right-hand side of the Salpeter equation (31). The 
result is 

d 3 q 



27 

x 



(2tt)3 ^(P-QMQ). 
[^Q + m) 1 °U xl °(p , Q-m)] ab 
(j)q — m 2 + ie) (p 2 Q — m 2 + ie^j 



(33) 



Now we substitute Eq. (22) into Eq. (33) and carry out the the p° integral on the right-hand 
side using residue theorem. If we close the contour in the lower half-plane, then we encounter 
two poles: 



p° = A/m 2 +p 2 - M mcson /2 - ie « \ (— - es] - ie, 

2 \m J 

p° = A/m 2 +p 2 + M mcson /2 - ie 2m - ie. 



(34a) 
(34b) 
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The first pole corresponds to the positive-energy pole in the quark propagator, and the 
second pole corresponds to the negative-energy pole in the anti-quark propagator. The 
contribution from the residue at the second pole, is suppressed as v 2 ~ p 2 /m 2 ~ |cb|/^ 
relative to the contribution from the residue at the first pole. Evaluating the latter, we 
obtain 

d 3 q 



^ b{p) w / [ v{p - qmq) 



J Q + m) 7 n A7 fa - m)} ab \ p0=Q 

x 7 2 \ ' ^) 

4m2 { e * ~ k ~ i€ ) 

where we have neglected corrections of relative order v 2 by evaluating the numerator at 
p° = 0. Making use of momentum-space Schrodinger equation, 

( eB _g^( P ) = y_g_y( P - g) ^), pel 

we find that 

*1(P) « -^Kh + mh°U(TlxU>h (h - m)] Vb $(p). (37) 
It is easy to see, by substituting this result into the right-hand side of the Salpeter equation 
(31) and repeating the preceding analysis, that this is a self-consistent solution of the Salpeter 
equation through terms of relative order v. 

Finally, returning to coordinate space, we can write the equal-time BS wave functions as 

(Vc(P)\QQ(x° = 0)\0) « - 1 (^ Q + m)(P + 2m) 1 5 {^-m)^(x), (38a) 



<JMP,e)|QQ(x° = 0)|0) « —L—y Q + m )y> + 2m)i(j> Q -m)i){x). (38b) 

As we have mentioned, p now has only spatial components. The interpretation of pP in 
coordinate space is — iV J . Using P -p = and P ■ e = 0, we can rewrite the equal-time wave 
functions as 

(r) c (P)\QQ(x° = 0)|0) « 7 2 ")|(7 5 + ^f , 7 5 + i ^^,f , ]7 5 )^), (39a) 

(J/m e)\QQ(x° = 0)|0) » -^±= (V + i [ 7 °, fd + ^ {A /} + ^/PW) 

x^(a;), (39b) 

where we have retained only terms up to those linear in p, and we have used 

{ 1 a ,a^ u } = -2e a ^ /3 7 /37 5 , (40a) 
j>ft+ftj> = 2ie a ^ Pa P^ lpl 5 . (40b) 
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By following a procedure similar to the one that we have just presented, one can also 
obtain the BS wave functions at equal light-cone time (x + = 0). In this case, one integrates 
the Fourier transform of the BS equation (30) over p~, rather than p°. The contribution at 
leading order in v still comes from the residue of the positive-energy pole of the heavy-quark 
propagator. Now, however, the mass-shell condition 

(P/2 + p)° = ^m 2 + p 2 (41a) 

or, equivalently, 

(P/2+p)' (P/2+p) + = m 2 +p 2 ± (41b) 

is used to eliminate p~ , rather than p° in the residue. Hence, one can obtain the equal-light- 
cone-time wave functions from the equal-time wave functions by using Eq. (41) to express 
p 3 in terms of p + . We work out the details of this transformation in Sec. IV. 

Comparing the definitions of the light-cone distribution amplitudes in Eqs. (18), (19), 
and (21) with the equal-time BS wave functions in Eqs. (39), we obtain the correspondences 



Pa <-> i>(r), (42a) 
P P <-> ip(r), (42b) 
x~P T <-> p~tp(r) ~ — p + ip(r), (42c) 



and 



V± <-> ip(r), (43a) 

V <-> 0, (43b) 

V T <-> ip(r), (43c) 

x~Va <-> p~~ip(r) fa — p + -?/>(r). (43d) 

Using Eqs. (20) and (43b), we obtain 

V L <-> ^(r). (43e) 

The meaning of these correspondences is that we take the Fourier transform of the quantities 
on the right-hand sides, use the mass-shell condition (41) to eliminate p 3 in favor of z — 
(P + /2+p + )/P + , and integrate over p± to form the light-cone distribution at zero transverse 
separation. We will give the details of these steps in Sec. IV. 
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IV. LIGHT-CONE DISTRIBUTION AMPLITUDES 



In this section, we compute the light-cone distribution amplitudes in terms of the 
potential- model wave functions, making use of the correspondences in Eqs. (42) and (43). 

Let us consider first the cases in which a light-cone distribution amplitude corresponds 
to ipij-), as for Pa, Pp, V±, Vt, and Vl- We begin by taking the Fourier transform of tp(r): 

^(p) = J d 3 xe tp - x 4>(r). (44) 

From the definitions of the light-cone distribution amplitudes in Eqs. (18), (19), and (21), 
we see that the light-cone momentum fraction z is related to the momentum variable p as 

n-P n-P ' v ; 

where n = (0 + , 1~, 0j_) is a light-like vector whose spatial component is parallel to the three- 
momentum P of the bound state. At the level of precision of this calculation, there is some 
ambiguity in the definition of z, in that we can discard terms of relative order v 2 . Hence, 
we can write in the quarkonium rest frame 

= P + /2 + p + = P + /2+p+ _ P + /2+p + = E Q +p 3 
Z P+ P° P° + 2p° 2E Q ' [ ' 

where we have made use of the mass-shell condition (41). The quark energy Eq is defined 

as 

E Q = y/\p\* + ml (47) 

We take the last expression in Eq. (46) as our definition of z. This definition is identical to 
that which is used in BC. It has the desirable properties that < z < 1 and that z <-> 1 — z 
under the interchange of the quark and the antiquark. In this definition, m c is the pole mass 
of the constituent heavy quark. Note that we distinguish m c from the parameter m that 
appears in the bound-state equation for our potential-model wave function. For numerical 
calculations in this paper we take m c = 1.4 GeV. 

According to the definitions of the light-cone distribution amplitudes in Eqs. (18), (19), 
and (21), the light-cone distribution amplitudes are given by the BS wave functions at 
zero transverse spatial separation, as well as zero light-cone-time separation. Therefore, in 
momentum space, we obtain the light-cone distribution amplitudes by integrating over the 
transverse momenta of the potential-model wave functions. Carrying out this integration 
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and making the change of variables to the light-cone momentum fraction z, we find that the 
light-cone distribution amplitude <f>(z) is given by 

where the prefactor is introduced in order to respect the conventional normalization condi- 
tion 

/ dz(j)(z) = l. (49) 
Jo 

Hence, we have from Eqs. (42) and (43), 

P A (z) = P P (z) = V ± (z) = V T (z) = (f>(z). (50) 

Now we can simplify the expression for the light-cone distribution amplitude. The vari- 
ables p 3 and |p| can be expressed in terms of z and p± = as 

z - - I 

P 3 = , n 2 . Jp± + m 2 c , (51a) 



The Jacobian in Eq. (48) is obtained from Eq. (51a) and Eq. (51b): 

dp 3 _ \Jp\ + m\ __ a/|/?| 2 + m 2 c 



(52) 



dz 4[z(l-z)f 2 " 2z(l-z) ■ 

The angular integral in Eq. (44) is simple to evaluate. Furthermore, for a fixed light-cone 
momentum fraction z, we can make use of Eq. (51b) to re-express the integral over the 
transverse momentum p± in Eq. (48) as an integral over the 



CO 

•2 



d 2 p ± = 8nz{l - z) / _d\p\ \p\. (53) 
Then, the expression for the light-cone distribution amplitude reduces to 

2 POO /"OO 

= — T77TT / _d\pW\p\ 2 + m c / drip(r)rsm(\p\r). (54) 

^lP\y) Jy/d^) Jo 

The lower bound of the |p| integral is a function of z: 

Vd(z) = m c \ Z 2 1 (55) 

y/z{l ~ Z) 
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which is the right-hand side of Eq. (51b) at p± = 0. The z dependence of (p(z) appears only 
in the lower bound \/d(z) of the \p\ integral in Eq. (54). The symmetry <f>(z) = 0(1 — z) is 
manifest. Furthermore, it is clear that the maximum value of (f)(z) occurs at z — 1/2 and 
the minimum value (zero) occurs at z — or z — 1. 

Now let us turn to the light-cone distribution amplitudes Pt and Va, whose correspon- 
dences with the equal-time BS wave functions are given in Eqs. (42c) and (43d), respectively. 
In each of these cases, the Fourier transform replaces the factor x~ on the left side of the 
correspondence with a factor d/dp + . Consequently, the solution of the resulting differential 
equation for Pt or Va contains an integration over p + . The expressions for Pt and Va that 
follow from these solutions are not normalizable because this integration and the factor p + 
on the right-hand side of each correspondence lead to an ultraviolet divergence. Such ul- 
traviolet divergences are characteristic of amplitudes involving corrections of higher order 
in v. They can be rendered finite by treating them within an effective field theory, such 
as NRQCD. Rather than introduce this complication into the present discussion, we simply 
drop such contributions to the cross section. That is, we drop the contributions that arise 
from the light-cone distribution amplitudes Va and Pt- These contributions are suppressed 
as v 2 , and, as we shall see, they have only about a 30% affect on the cross-section calcu- 
lation for the model light-cone distribution amplitudes that are used in BC. Hence, these 
contributions do not have a significant effect on the discussion of the order-of-magnitude 
discrepancy between the calculation in BC and the NRQCD factorization calculations. 



A. Nonrelativistic approximation 

In the nonrelativistic approximation, one neglects the relative momenta of the heavy 
quarks in comparison with m c . In the light-cone formalism, this amounts to taking p± = 0± 
and z — 1/2. In this approximation, the Jacobian in Eq. (52) reduces to 

dp 3 



dz 



2m r . 



approx 



The corresponding approximate version of Eq. (54) is 

8m c f 



^approx 



» = *(1 " *) 



drip(r) cos a/ d{z)r 



(56) 



(57) 



The approximate Jacobian (56) was used in BC. For a given potential-model wave function, 
it leads to a narrower light-cone distribution amplitude than does the exact Jacobian (52). 
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Since, in this paper, we are generally working at leading order in v, it is unnecessary (but 
not inconsistent) to include the corrections of higher order in v that are contained in the 
exact Jacobian. However, it is convenient to use the exact Jacobian because it preserves the 
normalizations of the light-cone distributions, and we have used it to derive all of the numer- 
ical results that are presented in this paper. As a check, we have also carried out numerical 
calculations using the approximate Jacobian (56), and those calculations are qualitatively 
consistent with the conclusions that we draw from the calculations that are based on the 
exact Jacobian. 



B. Expressions for the light-cone distribution amplitudes for efficient numerical 
evaluation 

The double integral in Eq. (54) is difficult to evaluate numerically because of the oscilla- 
tory nature of the inner integrand and because the integrand of the outer integral decreases 
only slowly as |p| — > oo. We can improve the numerical accuracy of the integrations by 
splitting the integral into two pieces: 



(f)(z) = <f> a (z) +4>b{z), 



(58) 



where 



<t>a{z) = 



fl^(O) JWd(z) 



d\p\ \p\ / drijj(r) r sin(|p|r), (59a) 
Jo 

— [ d\p\ ( a/|p| 2 + m 2 c — \p\] [ dril>(r)rsm(\p\r). (59b) 

*W) LfdH) v J Jo 



The integrand in <f) a (z) contains the leading behavior of the integrand in Eq. (54), but is a 
simpler form in which the \p\ integration can be carried out analytically. The integral 4>b{ x ) 
contains the remainder, whose integrand converges more rapidly as \p\ — > oo. We carry out 
the integration over \p\ in (f> a (z), treating the divergent integration as a Fourier transform 
whose result is defined in the distribution sense. The result is 



<Pa (Z) = l + 



^(0) Jo 



drip(r) < a/ d(z) cos a/ d{z)r 



sin 



z)r 



(60) 



We use Eqs. (59b) and (60) to evaluate <f>(z) numerically. 
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It is useful, in checking the accuracy of numerical-integration methods that we use, to 
know the values of the light-cone distribution amplitudes for special cases in which they can 
be computed analytically. The light-cone distribution amplitudes 0, <p a , and (f>b are computed 
for the case of a pure Coulomb potential in Appendix A. Their asymptotic behaviors at 
z = (or z = 1) are also given. In addition, it is simple to compute the value of a at 
z = 1/2. In this case d(z) = 0, and Eq. (60) reduces to 

0a(l/2) = 1. (61) 



C. Light-cone distribution amplitudes 

In Fig. 2, we show the light-cone distribution amplitude that results from substituting our 
potential-model wave function into Eqs. (59b) and (60) and carrying out the integrations 
numerically, taking m c = 1.4 GeV. For comparison, we also show the model light-cone 
distribution amplitude that is used in BC: 



<foc(z) = C{V )z{\ - Z 



Zl 



(62) 



I - Az{\ - z){\ - v 2 ) _ 

where v 2 = 0.3 and c(0.3) ps 9.62. As can be seen, the model light-cone distribution 
amplitude of BC is very similar in shape to the light-cone distribution amplitude that we have 
derived from a potential model. However, there is a significant difference in the functional 
forms in the tails at z = 0,1. The asymptotic behavior of the BC light-cone distribution 
amplitude at z = is (f>Bc( z ) ~ c(v 2 )z^ v \ In contrast, the potential-model light-cone 
distribution amplitude that we derive behaves at z = as (f>(z) ~ 87c-2 1 ^ 2 /(vrm c ) [Eq. (A8)]. 

In Fig. 3, we show the light-cone distribution amplitude </> a pprox(z) that arises from 
Eq. (57), in which the nonrelativistic approximation has been used for the light-cone-fraction 
Jacobian. As can be seen from the figure, this approximation leads to a narrower light-cone 
distribution than does the exact expression. 



V. PERTURB ATIVE CORRECTION CONTAINED IN THE LIGHT-CONE DIS- 
TRIBUTION AMPLITUDE 

As can be seen from Eq. (51b), when (f>(z) is used in computing the cross section, it 
yields contributions that arise from both large and small QQ relative momentum p. The 
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FIG. 2: The light-cone distribution amplitude 4>(z) that is derived from the potential model in this 
paper and the model light-cone distribution amplitude <f>Bc(z) that is used in BC. As is explained in 
the text, (j)(z) is computed from the potential-model wave function by carrying out the integrations 
in Eqs. (59b) and (60) numerically, taking m c = 1.4 GeV. 4>bc{z) is given in Eq. (62). 

large-p contributions are not treated correctly in a potential model, which is inherently non- 
relativistic. In the NRQCD factorization approach, this difficulty is dealt with by factoring 
the large-p contributions into short-distance coefficients and the small-p contributions into 
NRQCD matrix elements. The large-p contributions are then computed in perturbation 
theory, which is presumed to be valid, owing to asymptotic freedom. 

In the light-cone calculation, we also wish to identify, and ultimately remove, these large- 
st? contributions. In particular, we wish to remove that part of the large-p contribution that, 
in the NRQCD approach, is computed as an order-o; s contribution to the cross section, i.e., 
an order-o; s contribution to the short-distance coefficient. Of course, we could, in principle, 
remove contributions of order a 2 s and higher, as well. These contributions, presumably, are 
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FIG. 3: Comparison of the light-cone distribution amplitude </> a pprox(z) with 4>(z). </> a pprox(^) is 
computed from Eq. (57), in which the nonrelativistic approximation is used for the light-cone- 
fraction Jacobian. cf)(z) is computed from the exact expressions in Eqs. (59b) and (60). 

less important than the order-o: s contribution. 

We begin by writing the production amplitude for the heavy- quarkonium in a schematic 
form, which is depicted in Fig. 4: 



M = I ^W(p)^(p). 



(63) 



where ip(p) is the spatial part of the momentum-space wave function for the meson and 
H(p) is the short-distance production amplitude for the quark pair. Here we approximate 
the temporal part of the wave function with 2tc5(p°). That is, we neglect deviations of 
the quark energy from m, which are of relative order v 2 . In the rest frame of the meson, 
p is the momentum of the charm quark and — p is that for the antiquark. At this point, 
we could simply impose a cutoff on the integration over p in order to separate the large- 
and small-p contributions. However, we wish to maintain consistency with dimensional 
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FIG. 4: Diagrammatic representation of M, as given in Eq. (63). The circle labeled H represents 
the hard part of the production amplitude, and the oval represents the quarkonium wave function. 




k P 
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FIG. 5: Diagrammatic representation of A4 — Mo, as given in Eq. (65). The circle labeled H 
represents the hard part of the production amplitude, and the oval represents the quarkonium 
wave function. 

regularization, which is conventionally used to compute the short-distance coefficients in 
perturbation theory in the NRQCD method. Therefore, we use the method of regions [30] 
to compute the large and small-p contributions. For small p, we approximate H(p) by H(0). 
Then, we obtain the standard leading-order contribution to the cross section: 

d 3 p 



M< 



(27T)" 



f V(p)# (0) = V(0)i7(0). 



(64) 



We subtract A4q from A4 . This removes a scaleless power infrared divergence, which vanishes 
in dimensional regularization. The remainder is infrared finite, and it is dominated by large 
p. Therefore, we can use the Bethe-Salpeter equation with the Coulomb-exchange kernel 
to write ip(p) approximately in a form in which one Coulomb loop is exposed. The result, 
which is depicted in Fig. 5, is 

n l n 

I I 



M-M 



(27T) 



X 



—I 



(27T) 



r [H(p)-H 



(65) 
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where g s is the strong coupling, and T a is an SU(3) color matrix. Note that we use the 
nonrelativistic expressions for the Q and Q propagators because we wish to identify the 
contribution to the high-momentum tail of the light-cone distribution that arises in the 
potential model. In the potential model, the heavy quark and antiquark are treated in the 
leading nonrelativistic approximation. Taking the small-A; approximation for the integrand 
of the k integration, we obtain the order-o; s short- distance contribution to M: 

AM = / i^r H&n (+w) - I , - ~J , [Hip) - h(o)\. 

(66) 

We emphasize that AM is not arbitrary, but has a precise definition in NRQCD: 
AM/ip(0) is the order-o; s short- distance coefficient of the NRQCD operator tpx^ m a QCD- 
like theory in which the Q and Q are treated in the leading nonrelativistic approximations 
and the only interactions are Coulomb-gluon ladder exchanges. That theory is, of course, 
identical to the potential model. AM/?p(0) is equal to the diagrams in Fig. 5, with the wave 
function truncated and k set to zero. That equality is precisely the NRQCD matching condi- 
tion: The first diagram in Fig. 5 is the order-a,, correction to the on-shell hard-scattering in 
the full theory, and the second diagram in Fig. 5 is the order-o; s correction to the NRQCD 
operator times the leading-order short-distance coefficient. (We assume that the on-shell 
matrix element of iftx , which multiplies AAi/ip(0), is normalized to unity.) 

The integral over p° in Eq. (66) can be carried out by the contour method. Here, we 
assume that any singularities in H{p) are of order m away from the origin and, therefore 
correspond to contributions that are suppressed as p 2 /m 2 ~ v 2 relative to the leading 
contributions, which come from the poles at p° = ±p 2 /(2m). The result is then 

AM = Ana s C F m J [H(p) - H(0)} 

= 4na s C F m J { fj rip T ]+ H (P)' ( 67a ) 

where 

H(p) = H(p)\ p o =0 . (67b) 

Here, we have neglected deviations of p° from in H(p), with errors of relative order v 2 . 
The second equality in Eq. (67a), defines the + distribution with respect to p. We use the + 
distribution to write AM as a convolution of the correction to the momentum-space wave 
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function A^(p) with the hard-scattering amplitude H(p): 

d 3 p 



AM 



(27T) 



where 



Aip(p) = 4ira s C F ip(0)- 



m 



(68a) 



(68b) 



Therefore, we can remove contributions to the production amplitude that correspond to the 
standard order-o; s corrections to the NRQCD short-distance coefficient by modifying the 
momentum-space wave function according to ip(p) — > i>(p) — Aip(p). 

Similarly, we modify the light-cone distribution amplitude <j)(z) according to (f>(z) — > 
4>(z) — A(f)(z). The definition of A<f>(z) follows from the definition of <f>(z) in Eq. (48): 



A<P(z) 



1 f d 2 p ± dp 3 
^(0)7 (2tt) 3 ~lh 

4:7ima s CF 



f d 2 p ± 


dp 3 1 


1 (2tt)3 


dz p 4 



+z 



The + distribution with respect to z in the last line of Eq. (69) is defined by 



dz [f(z)] +z H(z) 



dzf(z)[H(z)-H(l/2)}. 



(69) 



(70) 



The evaluation of A<f>(z) can be carried out in analogy with the derivation of Eq. (54) from 
Eq. (48). The exact Jacobian in Eq. (52) is substituted into Eq. (69) and the p± integral 
can be written as an integral over as in Eq. (53). Then, the expression for A<f>(z) can 



be written as a single integral: 



A<P{z) 



y s c F J' 



d\p\ 



d(z) 



in-. 



jr 1 



(71) 



+z 



Using the identity 



/ 



dx 



Vx^TT l Vx 2 + l + l 
2x 2 + 4 og V^TT-i 



(72) 



we obtain an analytic expression for A<f>(z): 

2q> s Cf ( m \ I \/z(l — z) 



A<f>(z) = 



7Y 



m, 



4(^-|) 2 



1 

+ 4 



+z 



log 



1 + 2^(1 -z) 

1 - 2y/ Z {l ~ Z) 



(73) 



+z. 



We note that, because of the definition of A<p as a + distribution, A<p has the property 



I dzA(f)(z)=0. 
Jo 



(74) 
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That is, the effect of subtracting A0 from is to shift contributions from z near or 1 to z 
near 1/2 without changing the overall normalization of the light-cone distribution amplitude. 

As can be seen from Eq. (73), A0 diverges near z = 1/2 as \/{z — 1/2) 2 . When A0 
is multiplied by a function f(z) and integrated over z, the leading linear divergence at 
z — 1/2 is canceled by the factor f(z) — /(1/2) that appears in the integrand by virtue of 
the definition of the + distribution. The subleading logarithmic divergence cancels because 
the corresponding integrand is odd about the point z — 1/2. 

Since the Coulomb potential controls the behavior of the wave function at large momenta 
(short distances), we expect the asymptotic behavior of A(j>(z) near z = (or z = 1) to be 
identical to that of 4>{z). It is easy to see from Eq. (73), that the asymptotic behavior of 
A(j)(z) near z = is 

A<j>(z) ~ — v^, (75) 
nm c 

where we have made use of Eq. (17). (The asymptotic behavior near z — 1 can be obtained 
from Eq. (75) by replacing z with 1 — z.) As expected, the asymptotic behavior of A<fi(z) 
in Eq. (75) is identical to that of <p(z) in Eq. (A8a). 

In Figs. 6 and 7, we plot A<f>(z), along with <f>(z). Since A0 is a + distribution [Eq. (73)], 
it contains a large negative contribution at z — 1/2 that is not shown in Fig. 6. As can be 
seen from Figs. 6 and 7, A<f>(z) shows the expected approach to 4>(z) at z = and z — 1. In 
Fig. 7, we show this behavior in more detail near z — 0. The convergence of <f>(z) and A<f>(z) 
at the end points has a significant consequence in the evaluation of the cross section. As we 
have mentioned, we must subtract A(f>(z) from <f>(z) in order to remove contributions that 
are conventionally calculated as part of the order-o; s contributions to the NRQCD short- 
distance coefficient. As we shall see, the cancellation between <f>(z) and A(f)(z) at the end 
points results in great reduction in the cross section. 

VI. LIGHT-CONE CROSS SECTION FOR e+e~ -» J/if> + rj c AT B FACTORIES 

In this section, we use the light-cone distribution amplitude that we have derived from 
the potential model to compute the cross section for e + e~ — > J/ip + r\ c at the B factories. 
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FIG. 6: A<p(z) and 4>(z), plotted over the full range of z. A(p(z) is given by Eq. (73) and <f){z) is 
computed from the potential-model wave function by carrying out the integrations in Eqs. (59b) 
and (60) numerically. Note that A</> is actually a + distribution, as is described in the text, and, 
therefore, it contains a large negative contribution at z = 1/2 that is not shown in the figure. 

A. Light-cone cross section 

We begin by recording the expression for the light-cone cross section from BC. In order to 
compare with the results of BC, we have tried to make use of the same approximations that 
were employed in that paper. (In some cases, it was necessary to deduce the approximations 
that were taken in BC by exploring the effects of various approximations on the numerical 
results.) In this section, we make use of the notation of BC. In particular, we use Mq to 
represent the heavy-quark mass, and we distinguish its numerical value from those of the 
parameter m in the potential model and the charm-quark pole mass m c . 

The light-cone cross section is derived by standard techniques [24, 31] and is given in BC 
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FIG. 7: Acp(z) and (f>(z), plotted over a restricted range in z. A</>(z) is given by Eq. (73) and 4>(z) 
is computed from the potential-model wave function by carrying out the integrations in Eqs. (59b) 
and (60) numerically 



as 



a(e + e- - J^ + Vc ) = ( -j^V \F VP (s)f, 



fi V (76) 

O \ -'-'beam / 

where p is the three momentum of the J/ip in the e + e~ center-of-momentum (CM) frame 
and v^s = 2i?bcam is the e + e~ CM energy. The form factor Fyp(s) is given by 



Fvp(s) 



32tt 
~9~ 



% 4 



(77) 



where M = (3Mj/^ + M Vc )/4 is spin-averaged S— wave charmonium mass and fy and /p 
are the decay constants of the J ftp and the r) c , respectively. The quantities |p| and q^ are 
given by 



\P\ = V^eam " M* , 

ql = {E heam + \p\) 2 . 



(78a) 
(78b) 
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I in Eq. (77) is defined by 

pi pi 5 

I = / dxi / dy 1 a a {k 2 ) y2fi(x,y), (79) 
Jo Jo i=1 

where fc is the momentum of the virtual gluon. The labeling of the light-cone momentum 
fractions for the J/ip, x±, x 2 — 1 — Xi, and for the ?7 C , y 1 , y 2 = 1 — yi, can be read from Fig. 1 
of BC. The integrands fi(x,y) in Eq. (79) are given by 

1 (1 - 2yQ V ± (x) P A (y) 

h(x - !,) = 2^ (80d) 

/5(X ' y) " 8 ^ ~M2~ J d^y~) • (8 ° e) 

The function d(x, y) in Eq. (80) originates from the gluon propagator in Fig. 1 of BC and is 
defined by 

d(x,y) = |l = (*, + A) ( Sl + A) , (81a) 

Z^V. (81b) 

The functions s(x) and s(y) in Eq. (80) are from the charm-quark propagator in Fig. 1 of 
BC and are defined by 

six) = Xl + | , (82a) 
yiV2% 

Z a M 2 

s(y) = yi + (82b) 

xix 2 q 2 

One might suppose that these expressions should contain two factors of Z k m for each factor 
of Mq, as is the case in the corresponding expressions in Refs. [10, 11]. Here, we wish to 
compare directly with the results of BC, and so we use the expressions that are given in 
that paper. The inclusion of two factors of Z^ for each factor of Mq in Eq. (82) has a small 
numerical effect on the cross section, which we will describe in Sec. VI C. Following BC, we 
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also take 



E - ^ 

-^beam ^ ' 



\P 



M 2 



beam 



2-E'beam 

g « Vs- 2M 2 , 
/c 2 « xiyig 2 + 2Mj. 



(83a) 
(83b) 

(83c) 
(83d) 



The various renormalization factors that appear in Eqs. (80) and (82) are defined by 

(84a) 



Z p — 




- 3Cf/6o 



a s (Ml 
Z m {i? = k 2 ), 



Z° m = Z m {^ = CT 2 ), 

where bo = 25/3. The running MS mass in Eqs. (84) is given by 



(84b) 

(84c) 

(84d) 
(84e) 



M Q (/i 2 ) = Z m (/i 2 ) M Q . (85) 
Following BC, we take the strong coupling a s to be given approximately by 

2\ 47l7&o 

a ^ ] " log[^/(200MeV)T (86) 

In Eqs. (84), a 2 is the square of the four-momentum of the charm-quark propagator in Fig. 1 
of BC. Again, we follow BC and use the approximate expression 

1 



1 + 



(87) 



v0--y). 

The renormalization factors Zj in Eq. (84) were introduced in BC in order to account 
for the evolution of the heavy-quark mass and tensor and pseudoscalar currents from the 
scale Mq to the scales of the hard interactions in the production process. We note that such 
renormalization factors do not appear in the NRQCD calculations of the cross section. As 
we shall see, they have a large numerical effect on the cross section. 
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We use the following numerical values, which were given in BC: s = 112 GeV 2 , 
f v = f P = 400 MeV, M Q = M c = 1.2 GeV. We also take Mj^ = 3.096916 GeV and 
M Vc = 2.98040 GeV, which imply that M = 3.06779 GeV. We note that, in the light-cone 
formalism, fy is related to the leptonic decay width of the J/ip as T(J/ifj — > e + e~) = 
(167TQ; 2 /27)\fv\ 2 /Mjj^. (See footnote 9 of BC.) Therefore, the value of fy is independent 
of the shape of light-cone distribution and, hence, is unaffected by the treatment of the 
high-momentum tail of the light-cone distribution. The value of fp is fixed by making use 
of the fact that f P is equal to fy, up to corrections of relative order v 2 , which are neglected. 

B. Numerical results 

Now we compute numerical values of the light-cone cross section (76), using the light-cone 
distribution amplitudes that we derived from the Cornell potential model. 

As we have already mentioned, we drop contributions that arise from the light-cone 
distribution amplitude Va, as these are of higher order in v. That is, we set the term / 5 in 
Eq. (80) to zero. For the model light-cone distribution amplitudes of BC, this has the effect 
of reducing the cross section from the value a(e + e~ — > J/ip + r] c ) ps 33 fb reported in BC to 
(j(e + e~ — > J/ip + rjc) ~ 23.47 fb. Thus, we do not expect the absence of this contribution to 
affect the order-of-magnitude discrepancy between the light-cone and NRQCD factorization 
results that we wish to address here. 

The remaining light-cone distribution amplitudes are given in terms of the light-cone 
distribution amplitude that is derived from the potential-model wave function [Eq. (50)]. 
Hence, we take 



where <f>(z) is defined in Eq. (54), and we use the equivalent forms in Eqs. (59b) and (60) for 
the numerical evaluation. Substituting these light-cone distribution amplitudes into Eq. (80), 



row of Table II and are denoted by a . If we correct the light-cone distribution amplitudes 
by subtracting the piece A0 [Eq. (73)] that corresponds to the order-o; s correction to the 
NRQCD short-distance coefficient, then we obtain the results that are shown in the fifth 



V T (x) = V L (x) = V ± (x) = <f>( Xl ), 
P P (y) = P A (y) = 0(yi), 



(88a) 
(88b) 



we obtain the values for the cross section <r(e + e 



J/ijj + T) c ) that are shown in the fourth 
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TABLE II: Values for the cross section a(e + e~ — > J ftp + %) computed using the light-cone dis- 
tribution amplitudes that are derived from a potential-model wave function. As is explained in 
the text, do denotes the unsubtracted cross section, a denotes the cross section in which contribu- 
tions that correspond to order-a s correction to the NRQCD short-distance coefficient have been 
removed from the light-cone distribution amplitudes, and as denotes the cross section in which the 
light-cone distribution amplitudes have been replaced by 5 functions. The notation "BC" indicates 
that the values of the renormalization constants Z^ or a s are those that were used in BC. The value 
M c = 1.2 GeV was used in BC; the values M c — 1.4 GeV and a s — 0.21 were used in Ref. [4]. The 
notation "BC" indicates that, in Eq. (82), we have inserted two factors of Z^ t for each factor of 
Mq, as is described in the text below Eq. (82). 



z 


BC 


BC 


1 


1 


1 




BC 


BC 


BC 


0.21 


0.21 


M c (GeV) 


1.2 


1.2 


1.2 


1.2 


1.4 


a (fb) 


22.51 


25.03 


11.59 


7.58 


4.95 


a(fb) 


7.77 


8.19 


4.30 


3.29 


2.48 


<J5 (fb) 


5.39 


5.59 


3.04 


2.48 


1.95 



row of Table II and are denoted by a. Finally, if we replace the light-cone distribution 
amplitudes 4>(z) by S(z — 1/2), then we obtain the results that are shown in the sixth row 
Table II and are denoted by as ■ 

C. Discussion 

As can be seen from Table II, if we take the values of Z i: a s , and M c that were used in BC, 
then the value of the unsubtracted cross section a that is obtained by using the potential- 
model light-cone distribution amplitude does not differ substantially from the value 23.47 fb 
that is obtained by using the model light-cone distribution amplitudes of BC (with the 
contribution of Va set to zero). This is not surprising, since, as can be seen from Fig. 2, the 
potential-model light-cone distribution amplitude does not differ greatly in shape from the 
model light-cone distribution amplitudes of BC. 

However, once the light-cone distribution amplitudes are corrected by subtracting the part 
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A(f) that corresponds to the order-ai s correction to the NRQCD short-distance coefficient, 
then the cross section is reduced significantly — by a factor of approximately three if the 
values of Z iy a s , and M c are those that were used in BC. 

A further large reduction in the cross section occurs if we set the renormalization factors 
Zi equal to unity. These factors represent an attempt to resum large logarithms of the 
heavy-quark virtuality in the production process divided by M% . However, it is not clear to 
us that these factors, as employed in BC, correctly resum such logarithms. The factors Z p 
and Z t are inserted in an attempt to account for the evolution of the light-cone distributions 
from the scale M% to the scale of the heavy-quark virtuality. However, these renormalization 
factors account correctly only for the evolution of the first moments of those distributions. 
The factor Z m evolves the heavy-quark mass in the heavy-quark propagator. However, some 
of the factors of Mq in Eqs. (81) and (82) should not be evolved because they do not arise 
from the propagator mass, but, rather, from the approximation that the external heavy- 
quark lines are taken to be on shell with mass Mq. Specifically, some of the factors of Mq 
that appear in the quantities d(x,y), s(x), and s(y) [Eqs. (81) and (82)] arise in this way. 

After we have set the renormalization factors Zi equal to unity, the cross section is reduced 
to 4.3 fb, which is comparable to the values 3.78 ±1.26 fb (Ref. [4]) and 5.5 fb (Ref. [5]) that 
were obtained in the NRQCD factorization calculations in leading order in a s . For purposes 
of comparison with the results of Ref. [4] , we show in the last two columns of Table II the 
effects of taking the values of M c and a s that were used in that paper. It should be noted in 
making this comparison that the NLO value of the wave function at the origin that was used 
in Ref. [4] corresponds to /y m 489 MeV. Use of this value for f v (rather than 400 MeV) 
would enhance the cross section by about a factor 2.2. 

In BC, it is claimed that a large enhancement of the cross section occurs in the light- 
cone calculation because it takes into account the finite width of the quarkonium wave 
function. However, we see from Table II that once the A<j) part of the light-cone distribution 
amplitudes has been subtracted, the cross section does not differ greatly from that which 
is obtained by using 5-function light-cone distribution amplitudes. We conclude that any 
large enhancement that arises from the finite width of the quarkonium wave function must 
be coming from the region that corresponds to the order-a s corrections to the NRQCD 
short- distance coefficient, i.e. the region of large heavy-quark momentum (> m c ). Within 
our potential-model calculation, it is essential to remove this region, as it is not treated 
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reliably under the nonrelativistic approximations that go into that model. In comparing 
with NRQCD factorization calculations, it is necessary to identify the contribution from 
this region, even if one makes use of light-cone distribution amplitudes with the correct 
high-momentum behavior, because otherwise its contributions would be double counted in 
the NRQCD factorization corrections of order a s to the production cross section. We note 
that the order-o; s corrections to the production cross section are indeed large [7]. 

The absence of a large enhancement of a relative to as is somewhat surprising, since it is 
known that, in the NRQCD factorization approach, there are large corrections of order v 2 
(Ref. [4]) that account partially for the effects of the finite width of the quarkonium wave 
function. However, the order-f 2 corrections arise from several sources: the dependence of 
the amplitude on the QQ relative momentum and the difference between the charmonium 
mass and 2m c in the phase space and in the relativistic normalizations of the states. Of 
these corrections, only the first type is sensitive to the finite width of the quarkonium 
wave function. It is 3.87(t> 2 ) relative to the lowest-order cross section [4]. Here, (v 2 ) = 
(Pi) / (m 2 ((9i)), where the matrix elements are taken in the J ftp or r\ c state, and the NRQCD 
operators 0\ and V\ are defined in Ref. [6]. (Their matrix elements are proportional to 
the square of the wave function at the origin and the wave function at the origin times 
the second derivative of the wave function at the origin, respectively.) Taking {Vi) / (Oi) = 
0.50±0.09±0.15 GeV 2 (Ref. [19]) and m c = 1.4 GeV, we see that the order-w 2 corrections that 
arise from the dependence of the amplitude on the QQ relative momentum are approximately 
99% of the lowest-order cross section. The difference between a and as is only approximately 
29%. However, it should be remembered that the light-cone formalism omits contributions 
that arise from the transverse part of the QQ relative momentum. (Note, though, that the 
transverse-momentum contributions are suppressed by a factor v 2 M^ ieson / E 2 eiirn relative to 
the longitudinal- momentum contributions.) We also point out that, in dropping terms that 
are proportional to the light-cone distribution amplitude Va, we have omitted corrections 
of order v 2 that, in the case of the model light-cone distributions of BC, amount to about 
a 30% correction relative to the cross section computed without them. In order to be more 
quantitative about the effects of these omissions, we have calculated the order-f 2 correction 
that is associated with the finite widths of the light-cone distributions by expanding the 
light-cone cross section around z = 1/2. The result is a relative correction 2.1(t> 2 ) 54%. 
Presumably, the remaining small difference between this correction and (a — as) /a of about 



33 



54% — 29% = 25% arises from corrections of still higher orders in v and from the fact that 
we have computed S4> only at leading order in a s . 

VII. SUMMARY 

The discrepancy between theoretical predictions in the NRQCD factorization approach 
and experimental measurements of the rate for the process e + e" — ■> J/ip + r\ c at the B 
factories presents a significant challenge to our understanding of the physics of quarkonia 
and, perhaps, to our understanding of QCD. In Ref. [8] (BC), it was suggested that the 
solution to this problem may lie in using the light-cone formalism to take into account the 
effects of the relative momenta of the heavy quarks and antiquarks in the quarkonia. The 
cross section calculated in BC is about a factor 6-10 times larger than the cross sections 
sections calculated at leading order in the NRQCD factorization approach [4, 5]. Both the 
NRQCD factorization approach and the light-cone method are believed to approximate QCD 
in the limit in which the hard-scattering momentum transfer is much larger than the heavy- 
quark mass or Aq CD . Therefore, the apparent discrepancy between the two approaches 
presents a further, purely theoretical, puzzle. In this paper, we attempt to address this 
theoretical issue on several fronts. 

First, we use a potential model to construct light-cone quarkonium distribution ampli- 
tudes that are a good approximation to the true light-cone distribution amplitudes of QCD. If 
one knows the static heavy-quark potential sufficiently accurately, then the potential-model 
allows one to determine the quarkonium wave function up to corrections of relative order v 2 , 
where v is the heavy-quark (or antiquark) velocity in the quarkonium rest frame, (v 2 ~ 0.3 
for charmonium). We make use of the Cornell potential [17], which, for an appropriate choice 
of parameters, yields a good fit to the lattice data for the heavy-quark potential — thus giving 
us confidence in the connection of the potential model to QCD. We use the Bethe-Salpeter 
equation to relate the potential-model wave functions to the definitions of the light-cone 
distributions, which are written in terms of the Bethe-Salpeter wave function. Then, we use 
on-shell kinematics for the heavy quarks to relate momentum variables to light-cone vari- 
ables. The latter approximation is valid up to corrections of relative order v 2 . Finally, we 
integrate out the relative QQ transverse momentum in the Bethe-Salpeter wave functions to 
form the light-cone distribution amplitudes. The light-cone distribution amplitudes that we 
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derive agree reasonably well in shape with the model light-cone distribution amplitudes of 
BC. The calculated cross section changes very little if we replace the light-cone distribution 
amplitudes of BC with ours. Thus, we conclude that the distribution amplitudes themselves 
are not the source of the theoretical discrepancy. 

A part of the light-cone result for the production cross section arises from the region of 
large QQ relative momentum (> m c ). This region is not treated accurately in our inherently 
nonrelativistic potential model. Furthermore, it consists of contributions that are part of the 
order-a,, correction to the hard-scattering coefficient in the NRQCD factorization approach. 
Hence, it is necessary to identify these contributions in order to avoid double-counting 
the order-a,, correction in making comparisons between the light-cone approach and the 
NRQCD factorization approach. After we subtract these contributions from the light-cone 
calculation, the production cross section is reduced by about a factor of three. 

The light-cone cross section in BC contains renormalization factors Zj for the heavy-quark 
mass and the pseudo-scalar and tensor point-like vertices. If we set them to unity, a further 
reduction in the cross section of about a factor of two occurs. These renormalization factors 
have no counterpart in a conventional NRQCD factorization calculation. They represent an 
attempt to resum large logarithms of the heavy-quark virtuality in the production process 
divided by the square of the heavy-quark mass. However, as we have explained in Sec. VIC, 
we are not convinced that these renormalization factors, as employed in BC, resum such 
logarithms correctly. It is important to understand whether a correct resummation yields 
a significant enhancement to the cross section and, if so, to incorporate it into NRQCD 
factorization calculations, as well as into light-cone calculations. 

Once the factors Zj have been set equal to unity, the subtracted light-cone cross section 
is comparable in size with the NRQCD factorization cross sections, and it differs by only 
about 29% from the light-cone cross section that is obtained by replacing the light-cone 
distribution amplitudes with zero-width distribution amplitudes (5 functions). 5 

We conclude that there is no conflict between the light-cone and NRQCD factorization 

5 Even in the case of <5-function light-cone distribution amplitudes, there are small differences between the 
light-cone and NRQCD factorization cross sections. Some of these differences arise from different choices 
for the quarkonium wave functions at the origin, a s , and the charm-quark mass and from the inclusion 
of quantum-electrodynamic corrections in Ref. [4]. Other differences, presumably, are due to the fact 
that the light-cone and NRQCD factorization cross sections become equal only in the limit of infinite 
hard-scattering momentum transfer. 
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results. Further, it seems that enhancements to the cross section that arise from making 
use of a light-cone wave function of finite width are rather modest and are not sufficient 
to remove the discrepancy between theory and experiment. That discrepancy remains an 
important puzzle. At this point, the most promising avenue for its resolution seems to be 
the inclusion of corrections of higher order in a s (Ref.[7]) and v. 
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APPENDIX A: LIGHT-CONE DISTRIBUTION AMPLITUDES FOR THE 
PURE COULOMB POTENTIAL AND ASYMPTOTIC BEHAVIOR AT z = 0, 1 

In this Appendix, we compute analytically the light-cone distribution amplitudes 0, (f> a , 
and (f>b for the special case of a pure Coulomb potential. We also determine their asymptotic 
behaviors at z = (or z = 1). These analytic results are useful in checking the accuracy 
of the numerical-integration methods that we use in the more general case of the Cornell 
potential. 

We drop the linear potential in the radial equation (6). In this pure Coulomb case, as is 
well known, Eq. (6) can be solved analytically. For the ground state, the wave function u(p) 
and eigenvalue ( are 




(Al) 



Cio 



A2 
4 



(A2) 
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Applying the relations in Eqs. (3), (5a), (5b), and (7), we obtain 



R(r) = 2 7 ^e 



(A3a) 
(A3b) 




From Eq. (A3a), it follows that 



Coulomb 



^(0)e" 7cr 



(A4) 



where the wave function at the origin is given by 



(A5) 



We next use Eq. (A4) to calculate the light-cone distribution amplitude in the pure 
Coulomb case. Substituting Eq. (A4) into Eq. (59) and integrating over r, we obtain 



where we have used <p = (p a + (pb- Carrying out the integration over \p\ in Eqs. (A6), we have 



Note that Eq. (A7b) can also be obtained by substituting the wave function in Eq. (A4) 
into Eq. (60) and integrating over r. 

We can use our results in the pure Coulomb case to compute the asymptotic behavior 
of the light-cone distribution amplitudes at z — (or z — 1). At z — 0, the corresponding 
longitudinal momentum is infinite [Eq. (51a)], and only the short- distance behavior of if){r) 
is important. That short- distance behavior is governed by the Coulomb potential. 6 Taking 

6 By differentiating Eqs. (59a) and (59b) with respect to z and setting z = 0, one can verify that, at 
z = 0, only the large p behaviors of the integrands for 4> a and are important and, hence, that only the 
short-distance behavior of ip( r ) is important. This approach yields a simple expression for the asymptotic 
behaviors of the derivatives of the light-cone distribution amplitudes with respect to z, which can then 
be translated into the asymptotic behaviors of the distribution amplitude. 




(A6b) 



(A6a) 




(A7b) 



(A7a) 
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the limit z — > in Eqs. (A7a) and (A7b) and using 4>b = — <p a , we find the asymptotic 
behaviors 

4>{z) ~ a (z) ~ ^z 1 / 2 , (A8a) 



7rm 



c 



~ ^, 3 / 2 . (A81,) 



The asymptotic forms in Eqs. (A8) provide good tests of the convergence in the high- 
momentum region of the numerical integrations that we use to compute (f>(z) and (f> a (z). 
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